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We determine the spectral photon yield from a hot dense plasma irradiated by VUV-FEL light in 
a Thomson scattering experiment. The Thomson signal is compared to the emission background 
mainly caused by bremsstrahlung photons. We determine experimental conditions that allow for a 
signal-to-background ratio larger than unity. By derivation of the Thomson and the bremsstrahlung 
spectrum from linear response theory we present a consistent quantum statistical approach to both 
processes. This allows for a systematic treatment of medium and quantum effects such as dynamical 
screening and strong collisions. Results are presented for the threshold FEL-intensity as a function 
of density and temperature. We show that the account for quantum effects leads to larger thresholds 
as compared to previous work. 
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I. INTRODUCTION 



Thomson scattering is a well established technique for 
experimental investigation of plasma parameters. Ex- 
amples can be found in Refs. 0, S S, II S @- Ob- 
servables like particle density, temperature, composition, 
and degree of ionization can be spatially and tempo- 
rally resolved by analysis of the scattering spectrum Q. 
Until recently, coherent light sources have been avail- 
able only for the visible and near UV part of the elec- 
tromagnetic spectrum. Due to small critical density 
n cr it = uj 2 e m c /e 2 10 20 cm~ 3 of free charge carriers 
for optical probes, the applicability of Thomson scatter- 
ing using coherent sources has been limited to targets of 
relatively low density. 



Glcnzer et al. 8, 9] have shown and explored the possi- 
bility of x-ray Thomson scattering in solid density tar- 
gets using the Ti He-a line at 4.75 keV as probe light 
[lg j . A new alternative emerged with the development of 
VUV-free electron lasers (VUV-FEL), providing pulses 
of coherent radiation in the far (vacuum-) ultraviolet. 
At the moment, the VUV-FEL at DESY Hamburg op- 
erates at 32 nm wavelength ^] , corresponding to 38 eV 
photons. With this coherent light source, dense mat- 
ter up to solid densities of 10 23 cm~ 3 can be penetrated, 
see Refs. [l^,ll3j. Under these conditions, the Thomson 
spectrum permits the determination of electron temper- 
ature and density directly from the position and height 
of collective resonances, i.e. plasmons, showing up in the 
scattering signal . First experiments will be performed 
in the near future at the VUV-FEL facility at DESY at 



A = 32 nm FEL wavelength, while in future stages of the 
project, wavelengths from 13 nm (VUV-FEL) down to 
0.1 nm (X-FEL) will be available. 

Due to the large number of free charge carriers 
at the temperatures and densities considered, thermal 
bremsstrahlung emission, resulting from inelastic free- 
free scattering, contributes significantly to the emission 
background. Therefore, experimental conditions such as 
scattering angles, spectral properties of the probe and 
the detector have to be chosen as to obtain a maximum 
signal-to-background ratio. Background is to be under- 
stood as bremsstrahlung radiation, whereas signal corre- 
sponds to the photons having undergone Thomson scat- 
tering. 

So far, classical formulas for the bremsstrahlung emis- 
sion level going back to Kramers 44] have been used to 
determine threshold intensities of the external source to 
overcome the background due to bremsstrahlung [l5l . ll6j . 
The Kramers result, given below in Eq. I|19fl . is de- 
rived from the assumption of Keplerian trajectories of 
the emitting electron in the Coulomb field of an ion 
and integration over all initial velocities weighted with 
the Maxwellian velocity distribution function. Quan- 
tum features are only accounted for in a semiclassical 
way: The velocity integral extends over velocities v, ful- 
filling the condition mv 2 /2 > fkj, i.e. the kinetic en- 
ergy has to be larger than the photon energy. Fur- 
ther quantum properties, such as the finite photon mo- 
mentum as well as the quantum mechanical nature of 
the scattering process are not accounted for. By com- 
paring the Thomson signal strength at the laser wave- 
length A = 14.7 nm to the bremsstrahlung photon yield 
calculated from Kramers formula^ Baldis et al. find 
threshold intensities of 10 13 W/cm for typical values of 
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free electron density n e 
k B T = 100 eV 16j. 



10 cm 6 and temperature 
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In this work, we evaluate threshold conditions (intensi- 
ties) using improved expressions for the bremsstrahlung 
spectrum. As usual, corrections to Kramers formula for 
bremsstrahlung are described by the so-called Gaunt fac- 
tor • I n the simplest approach it is obtained by taking 
into account collisions between electrons and fixed ions 
in Born approximation. In dense plasmas, many-particle 
effects as dynamical screening become important. More- 
over, strong collisions have to be accounted for. We show 
in this paper how the Gaunt factor can be derived from 
linear response theory [Isj in a general way. Within this 
framework, modifications of the emission spectrum be- 
yond Born aproximation can be included in a systematic 
manner 0] as will be discussed later. 

We then apply our formulas to determine the threshold 
intensities for a broad range of experimental parameters 
(wavelength, spectral properties of detectors, and differ- 
ent materials), relevant for future experiments at DESY. 
Furthermore we compare Thomson and bremsstrahlung 
photon yield over a finite spectral range. Thereby, and 
by taking improved expressions for the bremsstrahlung 
cross section, we show, that even higher thresholds have 
to be reached in order to obtain a Thomson signal above 
the bremsstrahlung level at least near the plasmon res- 
onances. These peaks are much lower than the central 
peak being essentially an ion feature. 

The present work is organized as follows: In the first 
section we review the basic physics of Thomson scatter- 
ing and bremsstrahlung and how they can be expressed 
in terms of the dynamic structure factor and the dielec- 
tric function, respectively. Since these two quantities are 
related to each other via the fluctuation-dissipation the- 
orem |2(J, we are able to describe both processes on a 
common and consistent basis. 

We then compare the emission level due to 
bremsstrahlung to the Thomson signal whose strength 
is proportional to the flux of incoming photons, i.e. the 
power density of the external source. Thereby, we find 
expressions for the threshold power density as a function 
of particle density and temperature. In the last section 
we discuss our results for various sets of experimental 
parameters relevant for future experiments at the VUV- 
FEL. 



II. THOMSON SCATTERING AND 
BREMSSTRAHLUNG 

The central quantity of interest is the spectral power 
density dP/dV dAdfi, i.e. the rate of energy radiated 
per unit scattering volume dT^, wavelength dA, and solid 
angle dtt. The total spectral power density is the sum of 
the corresponding quantity for every radiative process in 
the plasma. In this work we focus on Thomson scattering 
and bremsstrahlung, i.e. 
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To unambiguously identify the Thomson signal, we re- 
quire that the Thomson power density is at least equal 
to the bremsstrahlung level, 



d 3 P T h 
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dVdXdfl ~ dVdXdfl 



(2) 



The Thomson spectrum is given by the intensity of the 
probe laser Jl , and the Thomson scattering cross section 
d 2 OTh/doi;dn. To account for the finite spetral band- 
width of the detector, one has to convolute each power 
spectrum with a detector function G(A). In practice, 
this is only relevant for the Thomson signal since the 
bremsstrahlung spectrum is slowly varying in the rele- 
vant frequency region. The Thomson power spectrum 
reads 



dvdxdn 



dfl duj 



(3) 

We have introduced the response function R(X), where 
the bar denotes the convolution with the detector func- 
tion. Note that we assume an optically thin plasma, thus 
radiation transport is neglected. Also, due to the short 
pulselength of the VUV-FEL (20-120 fs), we neglect the 
heating of the plasma due to the probe beam. 

The bremsstrahlung spectrum does solely depend on 
the plasma parameters density and temperature, so 
that, with a suitable formula for the bremsstrahlung 
spectrum, which we will abbreviate by the notation 
d 3 P hT /dVdXdtt = j(X) in the following, Eq. J3J defines 
a threshold intensity 
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(4) 



We will now briefly describe how expressions for the 
Thomson scattering and the bremsstrahlung spectrum 
can be obtained from a common starting point, i.e. the 
dielectric function of the plasma. 



A. Thomson scattering 

The cross section for Thomson scattering in a plasma 
can be given in terms of the dynamic structure factor 
(DSF), S(k,uj): 
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(5) 



see Ref. [21|,|22j for details. The variables (fe,w) are re- 
lated to the transferred linear momentum and energy, 
respectively, while fco = too/c and k\ = |feo + k\ are in- 
coming and outgoing linear momenta of the laser field, 
respectively. (dcr(0)/dO)xh is the Thomson scattering 
cross section for the isolated scattering event. It is given 
by the Klein-Nishina formula, its derivation can be found 
in standard textbooks of quantum electrodynamics, e.g. 
Ref. El. 
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In the nonrelativistic limit, the Thomson cross section 
for unpolarized light is given by 
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(6) 



Here, r c — e 2 / '47reom e c 2 is the classical electron radius, 6 
is the scattering angle between by k\ and k$. 

Equation JSJ shows that plasma collective properties 
can be measured in a Thomson scattering experiment. 
However, this requires an accurate theory of the dynamic 
structure factor. The dynamic structure factor S(k, to) is 
closely related to the longitudinal dielectric function by 
the fluctuation-dissipation theorem |24| 



S(k,u>) 



Ime, 1 (k,uj) 



(7) 



7re 2 n e 1 — exp (— hw/k^T) 

This relation can be utilized by applying an appropriate 
approximation to the dielectric function. Due to the mass 
ratio m e /mi <C 1, Thomson scattering on ions can be 
neglected. Therefore, we will approximate S{k,uj) by 
S ee (k,u>), the electronic DSF. 



B. Bremsstrahlung 

Radiative free- free transitions of electrons, known as 
bremsstrahlung, represent the main source of emission in 
a hot, fully ionized plasma. However, the bremsstrahlung 
process requires the presence of a scattering partner that 
carries the recoil momentum. Kinematically, emission by 
electrons which scatter on ions is dominant. 

For a thermally equilibrated plasma, emission, char- 
acterized by the spectral power density j(uj), and the 
absorption coefficient a{oS) as the relative attenuation 
of the intensity of electromagnetic waves propagating in 
the medium per unit length, are linked by Kirchhoff's 
law [2^ j{u>) = L(uj)a(uj), with the Planck distribution 
L(cj) = ftio 3 / [4n 3 c 2 (exp(Huj/kBT) - 1)]. Note that we 
consider the long wavelength limit, thus the wavevector 
k does not show up as an argument of the quantities j{uj) 
and a(u). The absorption coefficient can be obtained 
from the imaginary part of the transverse dielectric func- 
tion according to 



u> Im et(to) 
c n(u>) 



(8) 



where the index of refraction n(u>) is also linked to the 
transverse dielectric function by 



nH = -L(Rc et H + | et H|) 1/2 



(9) 



In the long wavelength limit, the longitudinal dielec- 
tric function, which appears in the fluctuation-dissipation 
theorem l(7|l. and the transverse dielectric function coin- 
cide [2(| . Thus, relations J7J and JSJ) enable us to treat 
both Thomson scattering and bremsstrahlung on a com- 
mon basis, namely by using an appropriate theory for the 
longitudinal dielectric function. 



C. Consistent approximations 

The aim of this work is to compare both radiation 
processes, Thomson scattering and bremsstrahlung, cal- 
culated in a consistent approximation. The comparison 
has to be carried out between the contributions of cither 
process in leading order of density. 

Bremsstrahlung occurs in second order of the coupling 
constant oqed = 1/137 as can be seen from the transi- 
tion amplitude wj? r expressed by Feynman diagrams, 



giving the Bethc-Hcitler cross section, see Ref. |23j. 
The transition amplitude involves a longitudinal field 
(Coulomb field) i.e. a scattering partner, say an ion of 
effective charge Z e g. Thus, bremsstrahlung is naturally 
of second order in density. A free electron does not emit 
bremsstrahlung, unless collisions take place. The chal- 
lenge is then to accuratly describe the scattering process 
itself. 

Born approximation, as given by Eq. (|10|l does not de- 
scribe the correct behaviour of the bremsstrahlung spec- 
trum in the case of collisions involving high transfer mo- 
menta, so-called strong scatterings. These can be in- 
cluded by ladder-summation of all one-photon exchange 
processes, which leads to the t-matrix. Thereby, the 
electron-ion interaction is treated accuratly in all orders. 
In the Coulomb limit, one finds Sommerfeld's expression 
for bremsstrahlung, see Eq. (12211 below, and Ref. [27| . 
On the other hand, due to the long range behaviour of 
the Coulomb potential, Born approximation as well as 
the Sommerfeld result suffer from infrared divergencies. 
Here, the account for screening leads to convergent re- 
sults. For details, we refer to Ref. 0|- I n our calcula- 
tions we will use Sommerfeld's formula Eq. I|22|l . Since 
we are interested in the bremsstrahlung spectrum at high 
frequencies, i.e. in the vicinity of the FEL frequency 
(fku — 38 eV), no screening effects are considered here. 
Note that the use of the t-matrix does not increase the or- 
der in density, instead, it gives the accurate Gaunt factor 
in leading order of density. Further many-particle effects 
like self-energy and vertex corrections lead to higher or- 
der contributions in density and, therefore, are not con- 
sidered here. For a consistent treatment of the self-energy 
of the scattering electron, given by multiple scattering on 
ions and its impact on the bremsstrahlung spectrum, see 
Ref. [13 . 

In the case of Thomson scattering, we adopt the same 
strategy: We calculate the contribution from lowest order 
in density. Unlike bremsstrahlung, already the first order 
in density gives a finite contribution. The amplitude for 
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Thomson scattering in second order of the interaction 



fcl, U)\ 

P,E P 
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does not contain any longitudinal field and therefore be- 
haves regularly over the whole spectral range. 

Thus, the DSF, which describes the Thomson scatter- 
ing in the medium, is taken in random phase approxima- 
tion (RPA) that accounts for the dynamical screening. 
This is all the more important, as the Thomson scatter- 
ing cross section is evaluated in the vicinity of the plasma 
frequency where dynamical screening is the dominant ef- 
fect |23|. Collisions are neglected, since they occur in the 
next to leading order in density. 

To conclude, a consistent comparison between Thom- 
son scattering and bremsstrahlung in media asks for the 
cross section of either process in leading order in density. 
This consideration leads us to the RPA for the dynamic 
structure factor, and a t-matrix ladder summation for 
the bremsstrahlung spectrum. 



III. LINEAR RESPONSE THEORY 

Thomson scattering and bremsstrahlung can be deter- 
mined by the longitudinal dielectric function e; (fc, uj). We 
will now briefly outline how the longitudinal dielectric 
function can be obtained in the framework of linear re- 
sponse theory. A discussion of its fundamental aspects 
can be found in |18|. for various applications to optical 
properties of plasmas, see |2(]|. In general, the longitu- 
dinal dielectric function e/(fc, u) is given in terms of the 
dielectric susceptibility x cc <(fc,w), where c and d label 
the different species in the plasma and further quantum 
numbers such as spin, 



e t (k,u) = 



1 



1 - E CC ' V cd {k)xcc'{k,u) 



(12) 



V cc ' (fc) is the unscreened Coulomb potential between par- 
ticles of species c and d, 



V cc ,(k) = 



q c qc> 



(13) 



Within linear response theory, the Kubo formula re- 
lates the susceptibility (response function) to the current- 
current correlation function 26] 



k 



uq c q d 



(14) 



where /3 = l/k&T is the inverse temperature, f^o is a 
normalization volume. The current operator for particles 
of species c is defined as 



p 



k/2,cp+k/2,c 



(15) 



One expresses the current-current correlation function 
by a force-force correlation function (J^ c ; J^. c /)uj+i V , us- 
ing the time derivative of the current Jk, c — i [H, Jk,c] /ft- 
The force-force correlation function is more suited for a 
perturbative treatment, see Ref. (2i| and Add. lAl where 
also definitions and useful properties of correlation func- 
tions are given. 



IV. CALCULATION OF THE 
BREMSSTRAHLUNG SPECTRUM 



Since the bremsstrahlung spectrum is evaluated at fre- 
quencies far above the plasma frequency, we can use the 
high frequency limit of Eq. (|12f) in order to derive the 
absorption coefficient. For details, we refer to App. 
In the high fre quen cy limit, the absorption coefficient (JHJ 



is given by 
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CCqLJ 2 



Re ( Jn Jt 



0,ei u 0,e/u+ir) 



ImG jj(ui+irf) . 



ccquj" 

Here we introduced the force-force Green function 
G jj(uj), defined in Eq. I|A6(I . which can conveniently be 
calculated using Feynman diagrams, see Fig. Q](a). G4 



G A 






(a) 



(b) 



(c) 



FIG. 1: Diagrammatic representation of G jj{uj^). (a) full ac- 
count of all medium effects by a four particle Green function, 
(b) Born approximation, (c) t-matrix approximation. 



denotes a four-particle Green function that contains all 
interactions between electrons and ions [l^. I* 1 the long 
wavelength limit considered here, e-e and i-i collisions 
vanish due to momentum conservation. We perform a 
sequence of approximations by selecting certain diagrams 
contributing to G4. 



A. Born approximation 

As the lowest order contribution with respect to the in- 
teraction potential, Fig.^(b) shows the Born approxima- 
tion. This approach reproduces the well-known Bethe- 
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Heitler formula in the nonrelativistic limit: 



16Z c 2 ff nin / e 2 \ / P 



3c 3 



V 47reo / V 27rrn c 
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(17) 
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(18) 



i?p is the free particle energy, for nonrelativistic particles 
Ep = h 2 p 2 /2m c holds. Details of the calculation are dis- 
cussed in 28]. In Eq. i|17|) the effective ion charge Z e g 
is used to account for the screening of the charge Ze of 
the nucleus due to inner shell electrons. In Ref. [lfij . Z e g 
is calculated as a function of temperature and density 
in the framework of Thomas-Fermi theory, see Ref. [23 
for details. Alternatively, Z c g can be obtained in the so- 
called chemical picture by solving Saha's equation, see 
Ref. |23. Especially for high temperatures, the chemical 
picture gives more reliable results, as the ionization equi- 
libria between different ionization stages are satisfied. We 
will therefore use the chemical picture in this work and 
calculate the mean ionization level with COMPTRA04 
|30|. The Thomas-Fermi model will only be applied for 
the purpose of comparison to results obtained by Baldis 
et al. in Ref. 0. 

The logarithm in Eq. I|17|) is the nonrelativistic limit of 
the Bethe- Heitler cross section HH H3 • The drawback of 



this result is the divergence of the bremsstrahlung spec- 
trum in the limit ui — > 0. However, physically this is of 
no importance, since for low frequencies the index of re- 
fraction is different from unity and modifies the spectrum 
significantly. This is known as dielectric suppression |33| . 

It is common use to write formulas for the emission 
and absorption due to (inverse) bremsstrahlung in terms 
of Kramers classical result 0] 



&Z 2 s n c ni f e 



3m c c 3 \47reo 



,2 \ 3 



67rm f 



1/2 



exp(— f3huj) , 
(19) 

multiplied with a correction factor gg(u)), called Gaunt- 
factor |17|. which takes into account medium effects as 
well as quantum corrections, i.e. 



j(u) = j K (u>) -gs(u) 



(20) 



With Eq. (|18f) . the Gaunt factor in Born approximation 
reads 

tjl (u) = — eKp(/3nw/2)K (/3tkj/2) . (21) 

7T 



B. Strong collisions 

The Born approximation assumes free particles as the 
in- and out states in the scattering amplitude. Taking the 
true scattering states, i.e Coulomb wavefunctions, leads 
to the so-called Sommerfeld formula 1271. 
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and 
I, 
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/ Ei + VEJ 



x 2 F 1 l + l- ir) f ,l + l-ir] i ;2l + 2; 



i(Vi+Vf) 



cff i w -^i/2 irq + i + %)rq + 1 + ^)1 
r (21 + 2) 



Ef = Ei — fku 



(24) 



77?^ = Z 2 ff Ry / 'Ei/ f is the Sommerfeld parameter with 

the Rydberg energy Ry = m c e 4 /2(A-Ke h) 2 ~ 13.6 eV, 
2Fi(a, b; c; d) is the hypergeometric function |34^ and 



r(x) is the Gamma function. As shown in Ref. |19| . 
Sommerfeld's expression is also obtained by a t-matrix 
ladder summation with a statically screened Debye po- 
tential (Fig.n](c)) if the limit of vanishing inverse screen- 
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ing length (Coulomb limit) is taken. It gives the correct 
Gaunt factor in the low density limit, which is considered 
here. 

Fig. |3 shows the dependence of the Gaunt factor on 
the photon energy for electron temperatures k^T — 
10 eV, 100 eV and 1000 eV. Over the large energy in- 
terval shown, the Gaunt factor in either calculation for 
a fixed temperature does not vary much and is of or- 
der unity, thus the widely used approximation to set 
§s(u>) « 1 is justified for estimates of the emission level 
as done in Refs. 0,0- More specific, §s(ui) = 1 is a 
good approximation for photon energies comparable to 
the temperature hiojk^T ~ 1 and for temperatures in 
the vicinity of the ionization energy {k^T ~ Z 2 Ry for 
hydrogenic systems). If one of these conditions is not 
met, one should use the full Sommerfeld expression 11' L'l) 
or appropriate approximations, see the detailed discus- 
sion in Ref. j3^. For the VUV-FEL experiments, both 
requirements are satisfied only roughly: The laser pro- 
vides photon energies of tuv ~ 40 eV and the optical 
pump laser will excite the plasma to temperatures of 
k-&T — 10 ... 50 eV, which is of almost the same order 
of magnitude as the first ionization energies of aluminum 
(6 eV, 19 eV, 28 eV), or hydrogen (13.6 eV) 36], materials 
presumably used as target in the experiments. Thus we 
take into account the Gaunt factor in t-matrix approxi- 
mation (Eq. lf'2'21 ) in our calculations. 
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FIG. 2: Gaunt factor gs(u>) in Born approximation Eq. 12 It 
and using Sommerfeld's formula Eq. 1221 . respectively. Re- 
sults are presented as a function of the photon energy hui for 
various temperatures. 

In Ref. 19j, a consistent treatment of the im- 
pact of dynamical screening and strong collisions on 
bremsstrahlung based on the Gould-DeWitt scheme [33 
is given. It is shown that the high frequency behaviour 
of the Gaunt factor is dominated by the t-matrix contri- 
bution, whereas dynamical screening can be neglected as 
long as the considered frequencies are large compared to 
the plasma frequency. 



V. DYNAMIC STRUCTURE FACTOR 

The total dynamic structure factor is defined as [36 



S cc (k,U>) 



1 



2irN 



dte iut (p c (k,t)p c (-k,0)) ; (25) 



(...) denotes the ensemble equilibrium average. Here, 
our discussion is focused on the electron structure fac- 
tor S ec (k,uj). The fluctuation-dissipation theorem (JJJ in 
connection with Eq. I|12fl allows to express the electronic 
part of the dynamical structure factor via the electronic 
dielectric susceptibility Xee(k,u>), 



S cc {k,uj) = 



Im Xcc(k,uj) 
7rn e 1 — exp (— hw/k^,T c ) 



(26) 



We will consider the classical limit, where the denomi- 
nator of Eq. (|26|l can be approximated by Hu /k^T c . In 
RPA, the following equation for the susceptibility holds 

M- 

X^ A (k,u;) = Xc(k,u;)5 cc , + x° c (k,u>)V%(k,u) X c ,(k,u) . 

(27) 

V sc (fc, ui) is the screened interaction potential which acts 
between particles of species c and c'. It satisfies the equa- 
tion 

V£(k,w) = V cc ,(k) + Y,Vcd(k)j&(k,w)V£(k, U ) . 
d 

(28) 

The free susceptibility Xc(k> w ) is obtained from Eq. 114fl . 
taking the current-current correlation function only in 
zcroth order with respect to the interaction, i.e. 



x ° c (k,u J ) = if3n Q - 

= 5 cd 



k' 2 



uq c q c 



I jz . jz \0 



f, 



p+k/2 fp-k/2 



E c p+k/2 -EC p _ k/2 -h{u J + l r 1 ) 



(29) 



fp is the momentum distribution function of species c. 
Again, a small but finite imaginary frequency irj, rj > 
has to be introduced, in order to fix the sign of the 
imaginary part and to obtain convergent results. In the 
case of a classical two-component plasma with different 
temperatures T c , we take the Maxwell distribution. In 
the limit h — > the resulting susceptibility reads 



XT\k,oj) = -n n c W(x c )/k B T c 



with the plasma dispersion function 
W{x c ) = l-2x c e- x - J e t2 dt-iV^a 



(30) 



(31) 



and the dimensionless variable x c — y/ ui 2 m c /2k 2 k B T c . 
Eqs. I|27|l and (|28|l can be solved algebraically for the 
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RPA-susceptibilities x^,? A (fc,a;) as shown in Add. IbI see 
also Ref. |3|. Evaluation of Eq. © with x* PA (fc,^), 
see Eqs. (GUJ), and (|B5)l . yields the electronic DSF in RPA, 
i.e. 



the position to compare the signal from Thomson scatter- 
ing to the bremsstrahlung background thereby treating 
both processes on a common basis. 



S ee (k,uj) 



1 + ZesaHTe/TjWixi) 



1 + a 2 W{x c ) + a 2 Z off (Te/T i )iy(a;i) 



Z, 



cir 



-a 2 W(x c ) 



x e exp (-xf) 



2 Xj exp (—x 2 ) 



1 + a 2 W(x c ) + a 2 Z eff (T e /7i)W(zO 



Here, the scattering parameter 



= n c /k = y/n e 2 /m c k B T e /k 



(32) 



(33) 



has been introduced. It separates the regime of collective 
Thomson scattering (a> 1), where collective resonances 
(plasmons) appear in the spectrum, and the regime of 
noncollective Thomson scattering (a <C 1), where the 
spectrum reflects the single-particle distribution function 
of electrons without further structures. Expression 1|32|) 
was also used by Baldis et al. in [l6| . 

Eq. lol't for the DSF contains two terms: The first 
term is basically the DSF of a free electron gas. How- 
ever, it also contains the ionic dispersion function, which 
accounts for the screening effect of ions. The second term 
gives the scattering signal from electrons which are close 
to ions and are therefore determined by the dynamics of 
these heavy particles. This term is dominant at small u) 
where W(x D ) can be approximated by its static limit, i.e. 
W(x c « 1) -> 1. For larger values of u, the first term 
dominates the spectrum, because the ion part is damped 
out more rapidly than the electron part. In particular, 
plasmon resonances, which appear in the spectrum due to 
a vanishing real part of the denominator, only survive in 
the first term, while they are damped out in the second. 
In so-called Salpeter approximation j4Qj the DSF (|32|1 is 
separated into two terms depending only on electronic 
and on ionic variables respectively, electron-ion correla- 
tions are neglected to a large extent. Therefore, the full 
expression (|32H should be used, as done throughout this 
work. 

Chihara 23] gives the DSF in terms of the local field 
correction factor, thereby including electron-ion colli- 
sions. In appropriate limits, his expression coincides with 
Eq. (J22J, see Ref. UJ. On the other hand, due to the 
equivalence of the local field correction and the collision 
frequency in the high frequency limit, Chihara's expres- 
sion leads to the same formula for the absorption co- 
efficient a(u) as given in Eq. JSJ 41] . Thus, Chihara 
presents an alternative approach to the question of pho- 
toabsorption and emission and Thomson scattering in 
plasmas starting from the DSF, whereas in this work the 
dielectric function is the central quantity. 

Having the thermal emission spectrum l|2(J|) and the dy- 
namic structure factor i|32|) at our disposal, we are now in 



VI. RESULTS 

We now compare power spectra for Thomson scatter- 
ing and bremsstrahlung emission. Figures |21 - El show 
results for different combinations of electron density and 
temperature for a hydrogenic plasma. We consider den- 
sities n e = 10 20 cm -3 and 10 22 cm -3 and temperatures 
of 10 eV and 50 eV as examples for possible experi- 
mental conditions with cryogenic targets. Furthermore, 
backscattering geometry (scattering angle 120°) is chosen 
and a VUV laser wavelength of A = 32 nm is assumed. 

To model the detector, we assume a detector function 
G(A) of gaussian shape, its width give n by the relative 
spectral bandwidth AA/A = Auj/ui 



givei 

m 



G(A) = —^c-* 2 / 2 " 2 



a = 0.425AA . 



(34) 



For the case of bremsstrahlung, we neglect the effect 
of the finite bandwidth of the detector, since the cor- 
responding spectrum is only slowly varying with fre- 
quency. The effective ion charge Z e g, calculated with 
COMPTRA04 is close to Z cfi = 1 in the case of 
hydrogen for the present values of electron density and 
temperature. 

In Fig. |3[ the black dashed curve represents the pure 
Thomson signal, i.e. no convolution with the detector 
function G(A) has been performed. For the present pa- 
rameters n c = 10 20 cm~ 3 , and /cbT = 10 eV, the Thom- 
son spectrum contains a very narrow ion peak, situated 
at the laser wavelength, and two satellites, which can be 
identified as electronic plasmon resonances. The central 
peak dominates these electronic resonances by a factor 
of 10, approximately. The solid black curve is obtained 
by convolution of the pure Thomson signal (dotted black 
curve) with the detector function. The detector resolu- 
tion is set to AA = 3.2 nm, which corresponds to 1% 
of the central wavelength. Due to the broad detector 
function, the central peak is lowered and broadened such 
that the plasmon resonances do not show up as sep- 
arate structures anymore. Instead, they only provide 
sidewings in the spectrum. In Fig. 21 (n c = 10 22 cm~ 3 , 
and kftT e = 10 eV), the electronic peaks are totally ab- 
sorbed in the central peak. In this case, the ion peak 
in the unconvoluted signal (not shown) is several or- 
ders of magnitude larger than the electronic resonances, 
which therefore do not contribute to the convolution in- 
tegral Fig. O shows spectra for n c — 10 20 cm~ 3 , 
and k B T e = 50 eV. At these parameters, one gets 
a = 0.56 < 1 for the scattering parameter, we are in the 
noncollective regime. Already the unconvoluted Thom- 
son signal is free of sharp electronic resonances, only two 
shoulders appear in the broad, noncollective signal. Fi- 
nally, in Fig.E0(n c = 10 22 cm" 3 , and k B T c = 50 eV), we 
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FIG. 3: Emission power spectra for Thomson scattering and 
bremsstrahlung. The black dashed curve represents the un- 
convoluted Thomson spectrum. Parameters: UbT c — 10 eV, 
n c — 10 20 cm -3 , laser wavelength A = 32 nm, scattering pa- 
rameter a — 1.25, Z c fi — 1 . 



FIG. 5: Emission power spectra for Thomson scattering and 
bremsstrahlung. Parameters: kuT c — 50 eV, n c — 10 20 cm -3 , 
laser wavelength A = 32 nm, scattering parameter a = 0.56, 
Z c s — 1 • 
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FIG. 4: Emission power spectra for Thomson scattering and 
bremsstrahlung. Parameters: feTL = 10 eV, n c = 10 22 cm -3 , 
laser wavelength A = 32 nm, scattering parameter a 
Z c s = 1 ■ 



12.5, 



FIG. 6: Emission power spectra for Thomson scattering and 
bremsstrahlung. Parameters: kuT c — 50 eV, n c — 10 22 cm -3 , 
laser wavelength A = 32 nm, scattering parameter a = 5.6, 
ZcS ~ 1 • 



have the same structure as in Fig. no plasmon peaks 
are visible due to the dominance of the ion feature. 

In Figures [3] - the intensity of the laser is cho- 
sen such that the central peak is situated clearly above 
the bremsstrahlung level. Three approximations for 
bremsstrahlung are shown, namely Kramers formula 
(|19fl . Born approximation i|18|) and Sommerfeld's formula 
l[22jl. Born approximation gives larger deviations from 
Kramers result than the Sommerfeld expression. This 
corresponds to the behaviour of the Gaunt factor shown 
in Fig. |21 The Sommerfeld result is always closer to unity 
than Born approximation. 

In the following, we evaluate the threshold intensity 
^thresh M defined in Eq. <gj) using Eqs. ©, ©, and (32J| 



for the Thomson power spectrum and Eq. (|20|l with the 
Gaunt factor in either Born approximation (|21|l or t- 
matrix approximation (|22|l for the bremsstrahlung power 
spectrum. 

Figures 171 - 1 1 01 show contour plots of the threshold in- 
tensity in part of the n c — T c -plane. In Fig.[71we compare 
Kramers result (dashed curves) to Born approximation 
(dotted curves) assuming aluminum as target material, 
scattering angle 9 = 20°, AA/A = 10~ 4 . In this case, the 
effective ion charge Z cr r is calculated using Thomas-Fermi 
theory j2?| as was also done by Baldis et al. in Ref. (lt| . 
who used the same set of parameters. Their results are 
reproduced by using Kramers formula for bremsstrahlung 
(dashed curve). 

For low temperatures, where the Gaunt factor in Born 



approximation is smaller than unity at the considered 
photon energy, cf. Fig. [21 higher densities are accessible 
as compared to Kramers result. For high temperatures, 
the opposite becomes true, g B (w) > 1 leads to lower ac- 
cessible densities. 

In figures |H1 - EH three approximations have been cal- 
culated for the bremsstrahlung level. Besides Kramers 
formula (dashed curves) we show results for Born approx- 
imation (dotted curves) and Sommerfeld's formula (solid 
curves). Comparing the Sommerfeld result to Born ap- 
proximation, it can be noted that Born approximation 
gives larger corrections, while Sommerfeld's theory leads 
to smaller deviations from Kramers result. This is one 
important result of this work: Taking into account quan- 
tum effects in a rigorous way via Sommerfeld's expres- 
sion for the Gaunt factor leads only to small corrections 
of threshold intensities, while Born approximation tends 
to overestimate these effects. This underlines the im- 
portance to go beyond Born approximation. For mod- 
erate and high temperatures, the Sommerfeld result lies 
systematically below the Kramers result, due to the in- 
creasing Gaunt factor at high temperatures. 

Furthermore, we investigated the influence of other ex- 
perimental parameters, namely the laser wavelength, the 
material, and the spectral bandwidth of the spectrome- 
ter. The latter parameter turns out to be of great impor- 
tance: By comparison of Fig. with Fig- EI one observes 
that notably higher densities are accessible in the case 
of small spectral bandwidth (AA/A = 1CP 4 in Fig. |!5J) 
than for the relatively large bandwidth (AA/A = 1CP 2 in 
Fig.©. 

Comparing Fig. [7| to Fig. [51 the influence of the 
Z-number of the material becomes apparent: Low 
Z-materials produce less bremsstrahlung than high 
Z-elements due to the Z 2 -proportionality of the 
bremsstrahlung cross section, cf. Eq. (|19|l . On the other 
hand, the Z-dependance of the dynamical structure fac- 
tor largely cancels out for the free electron part of the 
structure factor while it is roughly Z/(l + Z) 2 for the 
ionic part, cf. Eq. 132|) . 

Finally, important differences are noted upon chang- 
ing the laser wavelength from 32 nm (Fig. |5J) to 13 nm 
(Fig. ITrj)l especially at low temperatures. Both data 
sets have been calculated using the same spectral res- 
olution (AA/A = 10 -4 ) and material (H). Looking at 
Fig. E| (A = 13 nm) at low temperatures the thresh- 
old contours are almost parallel to the density axis. For 
higher temperatures, both wavelengths 32 nm and 13 nm 
give nearly equal threshold intensities, which do not de- 
pend on temperature. This can be understood from the 
bremsstrahlung spectrum, which becomes nearly inde- 
pendent of frequency and temperature at low hw/k^T. 

Since the information about temperature and density 
is stored in the position and height of the plasmon res- 
onances, we will now focus on experimental conditions 
to be met in order to separate the plasmon peak from 
bremsstrahlung background. Results for the threshold 
intensity are given in Tab. [U The electron density is set 
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FIG. 7: Threshold curves in the density-temperature plane for 
Al. Parameters: laser wavelength 14.7 nm, 6 = 20°, AA/A = 
1(T 4 , T c /Ti = 2. 




1000 



FIG. 8: Threshold curves in the density-temperature plane for 
H. Parameters: laser wavelength 32 nm, 6 = 120°, AA/A = 
10" 2 , Te/71 = 2. 



to n c = 10 21 cm -3 , for n c = 10 20 cm" 3 and 10 22 cm -3 see 
figures|21-|ni Two wavelengths are considered, A = 32 nm, 
the momentary VUV-FEL wavelength at DESY, and 
A = 13 nm, envisaged wavelength in the near future. 
Since the latter wavelength allows for very efficient x- 
ray optics to be applied, the spectral resolution has been 
reduced to AA/A = 10~ 4 , while 32 nm wavelength al- 
lows only for AA/A = 10 -2 . Finally the scattering an- 
gle has been chosen such that the plasmon peak is on 
the one hand well pronounced but on the other hand 
shifted far enough from the central peak as to be re- 
solved by the spectrometer. For detailed discussion of the 
Thomson scattering spectrum at various plasma param- 
eters, we refer to Refs. The rightmost column of 
Tab. [I] gives the Gaunt factor in Sommerfeld approxima- 
tion (Eq. (|22|) 'l. Note that in the case of the 32 nm wave- 
length, the threshold intensity is increased by 14 — 22% 
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FIG. 9: Threshold curves in the density-temperature plane for 
H. Parameters: laser wavelength 32 nm, 9 = 120°, AA/A 
Hi \ Te/Ti = 2. 
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FIG. 10: Threshold curves in the density-temperature plane 
for H. Parameters: laser wavelength 13.0 nm, 9 — 120°, 
AA/A = 10~ 4 , T c /Ti = 2. 



due to the Gaunt factor. 



VII. CONCLUSION 

In this work we have studied Thomson scattering and 
bremsstrahlung emission in warm dense matter in the 
context of plasma diagnostic experiments to be per- 
formed in near future. Out of these two competing, 
Thomson scattering is supposed to serve as a probe for 
plasma parameters, but also bremsstrahlung gives an im- 
portant contribution to the photon yield from highly 
ionized plasmas. Thus, experimental conditions have 
to be determined that allow for a maximum signal-to- 
background ratio. Here, we focused on the laser inten- 
sity. To this purpose, expressions for the Thomson sig- 
nal, which is given by the dynamical structure factor of 



TABLE I: Experimental parameters and resulting threshold 
intensities for different electron temperatures and laser wave- 
lengths. Fixed parameters: Z = 1, n e = 10 21 cm -3 . 



fc B T c [eV] 


A[nm] 


9 


AA/A 


A p i[nm] 


/thresh [W/cm 2 ] 


9% 


10 


32 


120 


10" 2 


33.1 


1.37 ■ 10 12 


1.14 


20 


32 


120 


10~ 2 


33.2 


4.92 ■ 10 12 


1.17 


50 


32 


120 


10 -2 


33.4 


9.29 ■ 10 12 


1.22 


10 


13 


120 


10" 4 


13.3 


2.95 • 10 9 


1.03 


20 


13 


45 


10" 4 


13.2 


1.55 ■ 10 11 


1.04 


50 


13 


45 


10" 4 


13.2 


3.49 • 10 12 


1.06 



the plasma as well as for the bremsstrahlung spectrum 
have been derived from a common starting point, namely 
linear response theory. This approach allows for a sys- 
tematic treatment of medium and quantum effects, such 
as dynamical screening and strong collisions. 

We have applied our formulas to determine thresh- 
old intensities for the external photon source (FEL) as 
function of density and temperature as well as their de- 
pendence on further parameters, namely the laser wave- 
length, the spectral resolution of detectors, and the mate- 
rial. In the discussion we focused on the bremsstrahlung 
spectrum. It was shown that Born approximation over- 
estimates the effect of collisions, i.e. it leads to rela- 
tively high Gaunt factors if compared to Sommerfeld's 
expression. Sommerfeld's formula (1221 gives the correct 
Gaunt factor in a weakly coupled plasma. At tempera- 
tures k-oT c > 20 eV threshold intensities calculated with 
the Gaunt factor in Sommerfeld approximation are larger 
than those obtained by Baldis et al. in Ref . [0| , who used 
Kramers expression, i.e. gs = 1. 

The Thomson signal was analysed within RPA which 
gives the contribution in lowest order of density. For 
densities investigated here, collisions are also relevant for 
the Thomson signal and lead to a pronounced change of 
the electronic part in S(k,X). This was investigated in 
Ref. However, the account for collisions does not 

alter the results for the threshold intensities performed 
in this work, since the Thomson scattering signal was 
evaluated at the laser wavelength, were the ionic feature 
of the DSF dominates. 

Finally, we have demonstrated, that Thomson scatter- 
ing can easily overcome the bremsstrahlung background, 
if laser intensities of 10 s W/ cm 2 — 10 13 W/ cm 2 are pro- 
vided. From this point of view, Thomson scattering ex- 
periments for plasma diagnostics using VUV-FEL radia- 
tion seem feasible. 
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po is the equilibrium statistical operator. The time de- 
pendence of observables is taken in the Heisenberg pic- 
ture, using the system Hamiltionian 



APPENDIX A: CORRELATION FUNCTIONS 
AND GREEN FUNCTIONS 



Correlation functions for two observables A, B are de- 
fined according to 



(A;B) 
(A; -B) w +i,j 



drTr[A(-^r)B t po] 



k-q,d a k,d a p,c 



pkq.cd 



(A2) 



o Performing integration by parts, the current-current cor- 

°2 relation function (|16fl can be expressed through a force- 

dte i(u+i v )t (A(t);B) , (Al) force correlation function as 



Luq 



I3fl k 2 ( j* ;J *)2 



-zo;(J fc z ;J|) + (j|;^) 



uj-\-ir) 



uj-\-ir) 



(J, 



(A3) 



with j§ te =_i[H,J§ ie ]/H. 

Also, it is convenient to introduce a generalized col- 
lision frequency in analogy to the Drude rela- 
tion [2Cj e(w) = 1 — Wpj/((jj 2 + iu}v(u))) where uj" 2 ^ = 

J2 C n ce^./(eom c ) is the squared plasma frequency. 
By comparison with Eq. ljT2*)l using (J^; J|) = e Wpj/ (3V, , 
we establish an expression for the collision frequency in 
terms of correlation functions 



v(u) — 7T- lim 



(J^Jk) 



klu+iri 



(A4) 



Further details can be found in Ref. |20| . Making use of 
Eq. JSJ), the absorption coefficient can be expressed as 



is given in the form of a force-force correlation function, 
cf. Ref. [2(j. Thus, the absorption coefficient is directly 
proportional to the real part of the force-force correlation 
function, which itself can be determined using perturba- 
tion theory. 

The relation to the thermodynamic Green function of 
two observables is given by 



(A-B) 



1 



7T/3 J-oo UJ UJ + IT] — UJ 



IuiGab(uj + iv) 



(A6) 



Using Dirac's identity one obtains Re (A; B)^ +iri 
Im Gab(u + iv)/ UJ f3 , which directly leads to Eq. I|16|) . 



APPENDIX B: THE RPA-SUSCEPTIBILITY 



V 



Re v(u>) 



c (ui 2 — 2ulmv{ui) + \v(uj)\ 2 )n(uj) 



(A5) 



In the high frequency limit lu 3> u p i, the index of re- 
fraction is unity, the imaginary part of the collision 
frequency tends to zero and the collision frequency is 
small compared to the frequency to. Then, we can 
consider the approximation a(uu) — uj^Re v{uj)/uj 2 c = 

(3floKe (Jq, </ ) w+ ^/ceow 2 , where the collision frequency 



Eq. H27fl is obtained from the exact relation for the 
response function Xcc'{k,uj), 

Xcc'{k,Lj) =Xc(k,u)5 cc ,+J2xc{k,u)V c s Z(k,L J )xdc>(k,uj) , 

d 

( B1 ) 

by truncation after the first iteration, i.e. Xdd — Xdc' m 
the second term of Eq. IjBll) . Insertion of the screened 
potential Eq. I|28|l yields the closed equation 
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V cc ,(k) + V r aI (fc)xS(fc,w)VS5(fc,w) 



(B2) 



For a two component plasma, i.e. an electron-ion 
plasma c = e, i, we obtain by matrix inversion 



,RPA 



Xc°(l-X?Mi) 



,RPA 



i - x° c v cc ~ X °Vn + xfx° e (ViiV ee - V ci V ic ) ' 

(B3) 

_ xXx° 

Xi ° 1 - X°eVee ~ + X?x2(Wee " V ei V ie ) ' 

(B4) 

X^; PA and Xn PA are obtained by interchanging indices i 
and e in Eq. (|B4J| and Eq. (|B3|t . respectively. 

The imaginary part of the electronic susceptibility 
in RPA, x^f A (fc,w) (Eq. (JB3J) for a hydrogen plasma 
{Vcc = V\\ = —V C i = —Vi C = V) now evaluates to 



This expression can also be derived from kinetic theory, 
i.e. the perturbative expansion of Vlasov's equation for 
a two component plasma 0]. It was used in to 
compare the Thomson signal to the emission background 
caused by thermal bremsstrahlung photons. Generalizing 
for different electron and ion temperatures T c and 7] and 
inserting the explicit expression for the response function 
in RPA iJSDJl, one obtains Eq. 



Im X l PA (k,co) 



l-V(k)xKk,co) 



l-V(k) [x° e (k,co)+xf(k,cv)} 
V(k) X ° e (k,io) 



l-V(k) [x°(k,co)+x?(k,tu)} 



lmx° e {k,u) 



Imx° (k,u>) 



(B5) 
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